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Abstract 

We present a view of the non-extensive thermodynamics based on general composition rules. A formal 
logarithm maps these rules to the addition, which can be used to generate stationary distributions by 
standard techniques. We review the most commonly used rules and as an application we discuss the 
Tsallis-Pareto distribution of transverse momenta of energetic hadrons, which emerge from relativistic 
heavy-ion collisions. 

1 Introduction 

In non-cxtcnsive systems those thermodynamic variables, which usualy scale with system size 
(with volume V and particle number N), violate somehow this scaling. Generally the entropy 
and energy of an A^-particle system, Sn and En, are composed from individual quantities. Si 
and El, by counting for interaction and correlation corrections. Considering pair interactions, 
mediated by the pair-potential v(r) (which depends on the relative coordinates r = ri — r2), 
and a particular form of the two-particle phase space occupation, pi2 = f{pi)fip2)g{'r) with the 
pair correlation function g{r), approaching one for uncorrelated pairs, the following relations can 
easily be derived: 

Sn/N = Si- n Jglng (fr, 

En/N = Ei+n Jgv d'^r, (1) 

while the volume and the particle number are related via the mean density: V = N/n — J df^r. 
In the above expression Si — 2s j g{r) d'^r with 's — ~ J f(p) \nf{p) d'^p and Ei = 2e J g{r) d'^r 
with e = J f{p)K{p) d'^p, where K{p) denotes the kinetic energy of a single particle. In principle 
the correlation g{r) can be obtained from the interaction v{r) in a stationary state (if exists), but 
this is a very complicated and difficult calculation for most of the known physical systems. Non- 
extensivity occurs, whenever the specific ratios, like En/N and Sn/N diverge in the TV — s- cx) 
limit. One considers this limit at constant mean number density, n. 

It is easy to construct examples for non-extensive energy at extensive entropy: in the crude 
approximation, when g(r) is either zero (up to a characteristic short range length) or one (towards 
infinity), the entropy correction is zero and hence the large- A'^ entropy is extensive. However, 
with a pair potential of the form v{r) ~ r~", whenever a < d with d being the spatial dimension 
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(d = 3 for isotropic systems), the correction to the energy becomes iV-dependent in the large- iV 
hmit Jj: 

En/N k Ei+ const. TVi-"/'^. (2) 

In systems with long range correlations (with scale independence, like in some random networks) 
g(r) may differ from one even at large distances. In such cases the entropy also may pick up a 
non-extensive contribution. Another example may be given by confinement: the correlation g{r) 
approaches zero in this case for confined pairs at large relative distances, but the integral of 5' In 5 
may have a non- vanishing contribution which ~ if of power-law type - may diverge with the total 
volume logarithmically. 

The composition of small sytems into a big one (extensivity) and the composition of two large 
systems (additivity) are related problems. Non-extensive systems are always non- additive by 
using the original definitions for energy and entropy. In a quite broad class of cases, however, 
another additive quantity may be constructed. The existence of such a quantity relies on special 
properties (specifically on the associativity) of the composition rule. The mapping of a non- 
additive quantity to an additive one, the formal logarithm, usually contains parameters, which 
describe the degree of non-extensivity. This is the basis of the construction of non-additive 
entropy (and energy) formulas. 

This way a basic problem occurs for any non-extensive thermodynamics: how do large sub- 
systems equilibrate, whose thermal state is described not only by a temperature, T, but also by 
a non-extensivity parameter, say q. In particular, for the Aczel-Daroczy-Chrvat-Tsallis entropy 
formula, how does equilibration occur between different (q, T) systems? Given two preheated 
systems, does a common stationary distribution occur, will it be a Tsallis-Pareto distribution, 
and is the temperature, defined by this equilibration process, universal (absolute)? We attacked 
such questions in the framework of a particular parton cascade model, using non-extensive energy 
composition rules in a Boltzmann equation type simulation [5] recently. 

We note that the often cited entropy formula 3J , 

^T = y^EK'-^^)' (3) 
using normalized probabilities ^^Wi = 1, follows the special composition rule 

Si2= Si+S2 + il-q)SiS2. (4) 

For factorizing probabilities, w^^^^ — w^^'''Wj^\ it can be mapped to an additive rule for Sr — 
Y{St), given as 

- - In (1 + aSr) ■ (5) 
a 

Here we used the parameter a = (1 — g). The result is the well-known Renyi entropy. 




which is additive but still contains the extra parameter q. For this additive, and hence extensive, 
entropy formula the question towards the two-parameter equilibration also holds [3] . 
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2 Abstract composition rules generalize 
non-extensivity 

The infinite repetition of an arbitrary pairwise, iterable composition rule is an associative rule [5]. 
It is a mathematical property that associative rules always possess a strict monotonic function, 
called here the formal logarithm, in terms of which they can be expressed 6^. We denote an 
abstract pairwise composition rule by the mapping {x,y) h(x,y). The associativity of such a 
rule is expressed by 

h{h{x,y),z) = h(x,h(y,z)) (7) 
for X, y and z being real quantities. The general solution of the associativity equation Q is given 

by 

h{x,y)^X-^{X{x)+X{y)) (8) 

with X{x) being a strict monotonic function, the formal logarithm. It maps the arbitrary com- 
position rule h(x,y) to the addition by taking the X-function of eq.([5]): 

X{h{x,y))^Xix)+X{y). (9) 

Due to this construction there are generalized analogs to classical extensive (and additive) quan- 
tities; they are formal logarithms. As a consequence stationary distributions, in particular by 
solving generalized Boltzmann equations [7] , are proportional to the Gibbs exponentials of the 
formal logarithm, 

fix) = (10) 

A general non-additive entropy formula can be derived based on the inverse of the formal loga- 
rithm (inverting the / expoX function): 

[ fX-\-\nf). (11) 



The rule leading to the g-exponential distribution is given by h{x,y) = x + y + axy with the 
parameter a proportional to q — 1. In this case one obtains the formal logarithm as being 
X{x) — i ln(l -|- ax). This formal logarithm leads to a stationary distribution with power-law 
tail as the function composition exp o X on the power —f3: 

f[E) = -e-'- = 1(1 + aE)-'"'' . (12) 

'Zi 'Zj 

The corresponding non-additive entropy formula is constructed as the expectation value of the 
inverse of this function, L^^ o In of 1//: 

f = -/(/i--/). (13) 

J a a J 

The Renyi entropy is the formal logarithm of the Aczel-Daroczy-Chrvat-Tsallis entropy belonging 
to this composition rule. 
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Further examples for non-additive rules can be easily given. The power-rule, h{x, y) — 
{x'' + y'^y^'', leads to a strecthed exponential, f{x) oc exp(— /3a;''), in the stationary state. Ka- 
niadakis [5] suggested a composition rule, h{x,y) = X\J\ + n'^y'^ 4- yVl + k^x"^ with the corre- 
sponding formal logaritm being the inverse sine hyperbolic function, X(a;) = ■i-Arsh(Ka;). The 

stationary distribution, /(x) — -g- (kx + Vl + K^a:^) ''^ , develops a power-law tail for large 
The corresponding entropy formula is the average of X^^oln of 1// over the allowed phase space: 

^ ■' ■. (14) 



Sk ^ ^ J ~ sinh(K In /) = J 
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Regarding kx — p/mc, the formal logarithm is proportional to the rapidity. This would imply a 
stationary distribution like ea;p(— 1?/|) with i] being the rapidity. 

Finally we note that the Tsallis rule h{x, y) ~ x + y + axy is particular, being the most general 
symmetric second order formula satisfying h{x, 0) ~ x. 



3 Equilibration of large subsystems 

Seeking for a canonical equilibrium state we have to maximize a total entropy given by a general 
composition rule, S{Ei, E2), at the same time satisfying a constraint which is in the general case 
also non-additive: h{Ei,E2) is constant. For the moment we neglect the dependence on further 
thermodynamical variables; usually the particle number N and the volume V is regarded to be 
proportional and extensive. 

In the traditional case both the entropy and the energy are combined additively: S{Ei, E2) — 
S{Ei) + S{E2) and h{Ei,E2) = Ei + E2- In the general case by using corresponding formal 
logarithms the quantities Y{S) and X{E) have to be considered as additive. Since for associative 
rules the formal logarithm is srict monotonic, the maximum of the total entropy is achieved where 
Y{S) has its extremum. The general canonical principle is therefore given by 

Y{S) - PX{E) = max. (15) 

The parameter (3 at this point is a Lagrange multiplier. Applying this for the equilibration of 
two large subsystems, and assuming that the entropy of each systems depends only on its own 
energy, one arrives at the equilibrium condition 

y'{S{Ei)) . Y'{S{E2)) 1 

Comparing this with the general canonical form eq. ()15|) we obtain that (3 — 1/T, and T is an 
absolute temperature in the clasiscal thermodynamical sense. Its relation to the entropy, however, 
has been generalized. In particular for an additive entropy, but non-additive energy composition 
rule, one arrives at 1/T — S'{E)/X'{E). The relation of this quantity to the logarithmic spectral 
slope, l/Tgiope = —d\nf/dE leads to a practical tool for the analysis of particle spectra in 
experiments. For the Pareto-Tsallis disitribution it is given by Tgiopc = T/X'{E) = T{1 + aE) = 
T + (q — 1)E. The naive effort to extract a temperature from energy spectra of particles, as 
it is a widespread praxis in relativistic heavy ion studies, only works if q — 1, i.e. for spectra 
exponential in the particle energy. Otherwise an energy dependent slope, and a curved spectrum 
in the logarithmic plot has to be interpreted. 
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4 Spectral temperatures in relativistic heavy ion collisions 



It is helpful to describe shortly, how a temperature can be conjectured from observations on par- 
ticle spectra produced in relativistic heavy ion collisions. The detected particles have relativistic 
velocities and different masses. 

One intriguing way is to look at the transverse momentum, px-, spectra around mid-rapidity 
(in the center of mass system for equal colliding heavy ions). The different identified hadrons, 
mostly pions, kaons, protons and antiprotons, have to show that their abundance in the momen- 
tum space depends on their energy; this phenomenon at zero rapidity is the so-called mT-scaling. 
The transverse mass is given as mr — \Jrn?- at strictly zero rapidity this is the total 

relativistic energy. 

The analysis is made a little more involved by the fact that the source emitting the detected 
hadrons is not at rest. The most prominent feature is a transverse flow, with relativistic velocity, 
vt (and a corresponding Lorcntz factor 7t = 1/ — ff^ in units where c—X). The relativistic 
energy of a particle in the frame of the emitting source cell is given by the Jiittner variable: 

E = Upp'' = jTiTT-T cosh(y — rf) — jtvtPt cos{(p — $). (17) 

Here the four- velocity of the source and the actual four-momentum of the particle are parametrized 
by rapidity and angle variables: 

= (7t cosh ry, 7t sinhry, 7TWT cos $, 7TWT sin $), 

— (mr cosh y,mT smh y,pT cos ip,pT sin (fi). (18) 

We consider a thermal model for the particle spectra; then the yield is supposed to depend on the 
Jiittner variable E given by eq. p7| . Assuming a general distribution f{E) exp{—X{E)/T), 
which is monotonic decreasing, one finds its maximum at the minimum of E. This variable is 
minimal at the rapidity ymin = ^-iid angle ipmin = ^, giving 

Emin = frmT - jtvtPt- (19) 

This Lorentz-boosted transverse energy reaches its minimum at the transverse momentum value 
PT,min = rujTVT j leading to uiT^min = nijT and £'„iin = rn. The expansion around this minimum 
in the p-r-distribution is an effective Gaussian: 

e-(-'")/-«expf-(^P^^^). (20) 

Such spectra are plotted in Fig[T]for a typical transverse flow of vt ~ 0.6 for massless and massive 
particles with masses of m = T and m — ST. The curves show typical shapes for light meson 
and heavy baryon spectra occuring in relativistic heavy ion collisions. 

In fact, according to experimental findings at RHIC the observed particle spectra have to be 
corrected for a transverse flow in order to reach m^-scaling. On the other hand the formula 
near the maximum, eq. ipUl) . may shed some light to the classical Einstein-Ott-Planck discussion 
about the temperature of relativistically moving bodies from an unexpected corner of modern 
experimental observations. Namely by using the Gaussian approximation both the particle mass 
and the effective spectral temperature gain a Lorentz factor, 77. 



5 



Effect of mass on Boltzmann pj spectra at v / c = 0.6 
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Figure 1: General shape of pt spectra for massive particles in the presence of a bulk transverse flow with 
radial velocity component vt = 0.6. The Gaussian approximation to the maximum reveals both a Lorentz- 
enhanced mass, m* = m7T, and a Lorentz-enhanced temperature parameter, T^t. 

5 Non-extensivity in quark matter and in hadron matter 

We conjecture that the power-law tails observed in hadronic spectra may stem from the non- 
extensivity of the preformed quark matter, which hadronizes rapidly. We make a connection 
between quark and hadron spectra by the quark coalescence model. A coalescence of two quarks 
(actually a quark and an antiquark) into a meson produces a yield proportional to the following 
quantity: 




(21) 



Here we integrate over the relative momentum of the quarks with a coalescence factor, C((7), for 
which a simple model has been utilized [S] . For common momenta much larger than the relative 
one |P| 3> I^tI (otherwise the quarks do not coalesce!) on gets 




(22) 



In particular light hadrons made from massless quarks follow the quark-scaling rule: 



/hadron (S) « /"(^/")- 



(23) 



As a consequence particular properties of the non-extensive thermal model between quark and 
hadron matter also scale: Tmesons = Tbaryons = ^quarks for the temperature, while ^mesons — 1 = 
((Jquarks ~ l)/2 for mcsons and gbaryons — 1 = (^quarks " l)/3 for baryons. Experimentally these 
relations are still to be checked. These predictions of the non-extensive phenomenology meet the 
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curves from pQCD calculations smoothly, with the following surmised properties of quark matter 
at RHIC: T = 140 . . . 180 MeV, q = 1.22, vt = 0.6. [10]. 

Summarizing we have shown that non-extensive behavior can be mapped to additive properties 
of a formal logarithm of the original quantity in the general case of associative composition 
rules. That such rules necessarily arise in the thermodynamical limit is demonstrated in Ref.[S]. 
Using formal logarithms all the classical concepts and techniques can be applied to describe 
thermal equilibrium or to generate distributions accordingly. The thermal equilibration of large 
subsystems are subject to straightforward generalizations of the familiar rules. As an exaple we 
discussed certain particle spectra arising in relativistic heavy ion collisions from a thermal and 
non-extensive quark matter. These results qualitatively agree with experimental findings. 
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